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(54) Procede et echographe pour la determination de la viscosite et du gradient de pression dans 
un vaisseau sanguin 



(57) Procede de determination de la viscosite et du 
gradient de pression dans un vaisseau sanguin, com- 
prenant I'acquisition de n>2 valeurs de vitessesdu sang 
en correspondance avec un meme nombre n de rayons 
du vaisseau sanguin, determines ie long d'un diametre 
situe a une position axiale donnee, formation d'un vec- 
teur de vitesses du sang au moyen de ces n valeurs de 
vitesses du sang, et evaluation de ladite viscosite et du 
gradient de pression a parti r d'une transformation dudit 
vecteur de vitesses du sang avec lormation d'une rela- 



tion lineaire liant directement un vecteur des debits (y) 
au vecteur des derivees de vitesses (h) factorise par la 
viscosite (u,) et au vecteur de gradient de pression (5), 
et evaluation simultanee des deux valeurs a determiner 
de viscosite (u.) et de gradient de pression (5) a partir de 
cette equation directe. Les vitesses sont favorablement 
determinees par correlation de signaux ultrasonores 
foumis par un echographe. 

Application : Imagerie medicale, echographie me- 
dical e 
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Description 

[0001] L' invention concerne un procede de determination de la viscosite et du gradient de pression dans un vaisseau 
sanguin. L'invention trouve son application dans le domaine de fabrication des appareils echographiques a uRrasons, 

s appeles echographies, pour I'aide au diagnostic medical. 

[0002] Dans le domaine medical, lorsqu'un patient a un accident vasculaire cerebral ou cardiaque, il est necessaire 
de contr6ler la viscosity du sang, car ce parametre joue un rdle important dans le diagnostic de ('accident vasculaire 
et dans la definition d'un traitement urterieur approprie, par exemple un traitement a base de medicaments fluidifiants. 
Actuellement, pour connaTtre la viscosite du sang, il est necessaire de realiser des mesures sur des echantillons san- 

io guins obtenus par des prelevements sur le patient. II en resulte des inconvenients importants. Le resultat des mesures 
n'est pas immediat, les prelevements eux-memes sont invasifs done potentiellement contaminants. Un autre inconve- 
nient est que les prelevements doivent etre repetes period iquement pour controler I'effet du traitement ce qui risque 
d'endommager les vaisseaux du patient. On ne connait pas actuellement d'alternative a cette methode de determination 
de la viscosity du sang. Par moyen invasif on entend un prelevement du sang qui necessite par exemple I' introduction 

15 d'une aiguille dans le vaisseau sanguin pour emplir une seringue avec un echantillon de sang. 

[0003] 1 1 est deja connu du brevet EP 0 430 374 (US 5 1 03 826, Bonnefous, Apr. 1 4, 1 992), un appareil echographique 
a ultrasons muni d'un systeme de traitement de signaux ultrasonores pour determiner des profils de la Vitesse du sang 
dans une artere. L'appareil echographique comprend un transducteur emettant et recevant des signaux ultrasonores 
d'une region d'artere d'un patient, et un etage transmetteur connecte a la sonde. Le systeme de traitement estime la 

20 vitesse du sang dans I'artere, en fonction du temps et de la profondeur d' exploration, par un procede de correlation 
de signaux successifs renvoyes par le transducteur. 

[0004] Un but de l'invention est de fournir un precede apte a etre mis en oeuvre, sans utiliser de moyens invasifs, 
pour permettre la determination de la viscosite du sang s'ecoulant dans une artere sous Taction de la pulsation car- 
diaque qui cr6e un gradient de pression, cette operation etant effectue avec une precision egale ou superieure a celle 
25 qui etait atteinte par la methode connue incluant I'analyse de prelevements. 
[0005] Ce but est atteint par un precede selon la revendication 1 . 

[0006] Ce procede presente des avantages importants du fait qu'il est aussi precis que le procede par prelevements, 
qu'il est plus rapide et plus simple a mettre en oeuvre. et qu'il peut etre repete un grand nombre de fois sans endom- 
mager les vaisseaux sanguins du patient, puisqu'il n'est pas invasif. Un autre avantage important reside dans le fait 
30 qu'aucun echantillon de sang n'etant preleve ou traite, tout probleme de contamination du patient ou du personnel de 
laboratoire est evite. 

[0007] Un appareil echographique a ultrasons selon la revendication 10 est muni de moyens pour mettre en oeuvre 
ce procede. 

[0008] Le procede est deer it ci-apres en detail, en reference aux figures schematiques annexe es dont : 

35 

la FIG. 1 qui represente sous forme de blocs fonctionnels, un diagramme de deroulement des etapes du procede 
de determination de la viscosite du sang dans une artere et du gradient de pression du a la pulsation cardiaque ; 
la FIG.2A qui represente un axe d'exploration echographique (Z) vis-a-vis de I'axe longitudinal (X) d'un vaisseau 
sanguin ; la FIG.2B qui illustre le decalage temporel entre deux tirs echographiques et la FIG.2C qui illustre une 
40 fenetre de resolution; 

les FIGS.3A et 3B qui illustrent des premier et second signaux temporels et la FIG.3C qui illustre un signal de 
correlation ; 

la FIG. 4 qui illustre un signal de correlation en 1 bit; 

la FIG.5 qui illustre un profit de vitesse le long d'un axe de tir echographique Z ; 
45 la FIG.6 qui represente sous forme de blocs fonctionnels un diagramme de deroulement du procede de transfor- 

mation d'un vecteur de vitesses pour fournir les valeurs de viscosite et gradient de pression dans le vaisseau 
sanguin. 

[0009] L'invention concerne un appareil d'echographie ultrasonore muni d'un systeme de traitement des signaux 
50 ultrasonores pour la determination des parametres physiologiques du sang s'ecoulant dans une artere examinee in 
vivo dans une region tissulaire, et en particulier pour la determination de la viscosite du sang de maniere non invasive 
et en temps reel, et pour la determination du gradient de pression sanguine dans cette artere. 
[0010] En reference a la FIG.1 et aux FIGS.2A, B, C, cet appareil comprend au moins un transducteur ultrasonore 
100 associe a un etage d'emission 200 et a un etage de traitement numerique 300. Le transducteur 100 est applique 
55 h la region tissulaire 11 pour I'examen de I'artere 10 ayant des parois 12, 13. Le faisceau £mis par le transducteur 
ultrasonore est de preference large bande. 

[0011] L'6tage d'emission 200 comprend un sequenceur compose d'un oscillateur et d'un diviseur de frequence qui 
commande, avec une periode de recurrence T, un generateur dont les signaux electriques d'excitation sont envoyes 
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vers le transducteur 100 qui les convertit en signaux temporels ultrasonores, qui sont emis avec la periode de recur- 
rence T, selon un axe de tir Z. Dans ce systeme d'emission, chaque signal temporel ultrasonore a un rang d'emission 
dans le temps note k, k+1 ,..,de periode T. L'etage d'emission 200 relie au transducteur 100 assure ainsi la formation 
d'un faisceau 50 ultrasonore pour I'exploration de la region tissulaire 11 incluant I'artere 10. Dans cet examen, on 

5 suppose que recoupment sanguin se fait selon I'axe longitudinal X de I'artere 10, et que I'axe de tir Z du faisceau 
ultrasonore a. partir du transducteur 100 fait un angle 0 avec I'axe X. Les profondeurs d'exploration sont discretisees 
sur I'axe de tir Z. en des points notes Rq, R v ... R n rep6res k partir d'un point 0, de coordonnee spatiale sur Z notee 
r c> situe a I'intersection des axes X et Z a I'interieur de I'artere, jusqu'a un point R N de coordonnee spatiale sur Z notee 
r N et situe sur une paroi par exemple 12 de I'artere 10. 

10 [0012] En referenced la FIG.1, un etage separateur 101 est insert ent re les etages 200, 300 et le transducteur 1 00 
pour eviter I'aveuglement des etages 200 et 300 par les signaux d'emission. L'etage de traitement 300 comprend en 
sortie de l'etage separateur 101 un amplrficateur haute frequence incluant une compensatbn de gain en lonction de 
la profondeur d'exploration dans la region tissulaire examinee. L'etage de traitement 300 comprend aussi un eliminateur 
301 d'echos fixes qui foumit, a partir d'un signal S k (t,r) provenant de l'etage separateur 101 , un signal echographique 

15 temporel d k (r,t) debarrasse des composantes fixes dues notamment aux reflexions speculates sur les parois 12, 13 
de I'artere 10 examinee. 

[0013] Le signal d k (r,t) issu de I'eliminateur d'echos fixes 301 est traite de preference selon une methode de corre- 
lation decrite dans le brevet EP 0 225 667 (US 4 803 390 : Bonnefous). Cette methode est brievement rappelee ci-apres. 
[0014] Selon cette methode, le systeme de traitement 300 determine la vitesse d'ecoulement de structures du sang : 

20 par exemple des amas de globules rouges, appelees cibles, notees TG (en anglais : Target), a I'instant t ou ces cibles 
coincident avec I'axe de tir Z et en des points discretises sur I'axe de tir Z appartenant a I'interieur de I'artere. li est 
suffisant d'estimer les vitesses entre le centre 0 et le point R N sur la paroi 12 car les vitesses entre le centre 0 et le 
point R' N k I'intersection de I'axe Z et de I'autre paroi 13 de I'artere sont supposdes symetriques. Les profondeurs 
d'exploration le long de I'axe Z sont donnees par des coordonnees spatiales r. Le temps ou un echo est renvoye par 

25 une cible situee a une profondeur r est note t. En particulier, le systeme de traitement 300 traite les signaux echogra- 
phiques temporels de rang k notes d k (r,t) formes des intensites des 6chos renvoy6s par les cibles en chaque point 
discretise de I'axe Z de profondeur r 0 k r N dans I'artere 10, et fournit la vitesse d'ecoulement du sang en ces points. 
Les FIGs. 3A et 38 represented de tels signaux temporels echographiques d k (r,t) et d k+1 (r.t) avec en ordonnee I'in- 
tensite du signal echographique en unite arbitrage et en abscisse des coordonnees temporelles t en microsecondes 

30 qui correspondent aux coordonnees spatiales r le long de I'axe de tir Z. 

[0015] Les FIGs.2A et 2B illustrent que I'appareil echographique n'est pas sensible aux vitesses autres que celles 
qui se rapportent k i'axe de tir Z. Done, on ne considere ni les distances qui sont en realite parcourues par les cibles 
TG selon I'axe X, ni les vitesses de ces cibles selon cet axe X, mais seulement leurs projections sur I'axe Z. La vitesse 
d'une cible est done mesuree par sa composante V z (r,t) le long de I'axe Z en fonction du temps t et de la profondeur 

35 d'exploration r. Sur la FIG.2C, on a note r k et r k+1 les profondeurs d'une cible TG renvoyant des echos dans les signaux 
de rang k et k+1 respectivement. La composante selon I'axe longitudinal X se deduit de V 2 (r,t) par la relation : 



V x (r f t) = V z (r,t)/cos9 (1) 

40 

[0016] Le systeme de traitement 300 met en oeuvre une technique de correlation de deux signaux temporels echo- 
graphiques successifs notes de rang k et k+1. La technique de correlation demande, en reference k la FIG.2C : qu'un 
fort pourcentage de cibles present dans le faisceau ultrasonore lors d'un premier tir de rang k reste present dans le 
faisceau echographique lors d'un tir ulterieur de rang k+1 , de maniere a ce que la correlation entre les deux signaux 

45 temporels successifs sort suffisamment forte. C'est pourquoi, en reference k la FIG.2B, la correlation est effectuee 
dans une certaine fenetre de resolution notee RG de dimension w (en anglais: Range Gate) choisie pour que cette 
demande soit facilement satisfaite. Par exemple, si le sang a une vitesse de flux de I'ordre de 1 m/s, si la periode T 
entre deux excitations est egale k 100 u.s ce qui correspond a une frequence d'excitation de I'ordre de 10 kHz, le 
deplacement des cibles entre deux tirs successifs est de I'ordre de 0,1 mm. Ce deplacement est substantiellement 

50 petit devant les dimensions d'une cellule de resolution RG avantageusement choisies de I'ordre de 1 mm selon I'axe 
X et 0 : 3 mm selon I'axe Z. Dans ce cas, environ 20 % de cibles quittent la cellule de resolution pendant la periode T 
environ 20 % y entrent et un fort pourcentage de cibles reste present. 

[0017] L'idee basique du traitement mis en oeuvre par le systeme 300, pour la determination de la vitesse en un 
point discretise de I'axe Z, est qu'un signal temporel ultrasonore c^ +1 (r,t) retrodiffuse par une cible en mouvement, 
55 situee k I'interieur de la fenetre de resolution RG, qui a ete emis par un tir echographique selon Z, de rang k+1 , au 
temps t, est la replique d'un signal temporel anterieur emis par un tir echographique de rang k selon le m§me axe Z 
au temps t -T, retrodiffuse par la meme cible au temps t - T t, ou t est le temps supplemental appele decalage 
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temporel (en anglais : Time Shift) mis par le signal temporal ulterieur pour parcourir le trajet sonde-cible-sonde, du fait 
du deplacement de la cible durant I'intervalle T. 

[0018] Ainsi, au premier ordre, deux signaux temporels echograph iques success if s ret rodiff uses par la meme cible 
sont lies par la relation : 

d k+1 (t) = d k (t-T-x) (2) 



10 



[0019] D'ou il resulte que le decalage temporel x est donn£ par : 

t = 2V 2 (r,t)/c (3) 



ou c est la Vitesse du son dans le vide et par extension dans les tissus. La formule (3) permet de determiner les 
composantes V 2 (r,t) des vitesses des cibles en chaque point discretise d'ordonnees spatiales r en faisant une approxi- 
mation basee sur le fait que la Vitesse axiale V 2 (r,t) est tres petite devant la vitesse du son c, et que le deplacement 
d'une cible est tres petit devant ia distance du transducteur a la cible durant une p6riode T. 

[0020] Par consequent, la vitesse du flux peut etre mesur§e par un appareil echographique ayant un systeme de 
traitement 300 munis de moyens pour mesurer ie decalage temporel x induit par les cibles entre deux tirs echographi- 
ques. 

[0021] Au premier ordre, la fonction d'intercorrelation entre les deux signaux temporels successifs est de la lorme : 



<Wto,u) = £*"* d t(0 d k ^(t + u)dt (4) 

ou le temps t 0 est lie a la profondeur d'exploration r par tg= 2r/c, w 0 est la largeur de la fenetre de resolution sur I'axe 
du temps, et u est un parametre de decalage temporel de la fonction d'intercorrelation (en anglais : cross-correlation). 
[0022] En introduisant I'equation (2) dans I'equation (4), il vient : " 

C k>k+1 (t o .u) = C kk (t 0 ,u-x) (5) 

[0023] La fonction C kM (t c ,u - t) est une fonction ^autocorrelation, et, ce fait est maximale pour u = x . Dans une 
hypothese d'une situation idSale de vitesse axiale uniforme, le maximum de la fonction de correlation correspond au 
decalage de temps x et par consequent permet une determination de la vitesse du sang pour le tir echographique au 
rang k selon la relation : 



40 V z (r,t) = (c/2) (x/T) (6) 

En reference a la FIG.1, le circuit 302 d'intercorrelation du systeme de traitement 300 determine les vitesses V 2 (r,t) 
en cherchant pour quel parametre u la fonction C|< k+1 (t^u) est maximale. A cet effet, la fonction d'intercorrelation est 
45 echantillonnee a un pas d'echantillonnage At, de facon a obtenir 2m+1 valeurs de la fonction de correlation, m etant 
un scalaire. La valeur maximale de ces 2m+1 valeurs correspondant a u = u 0 permet de mesurer x en utilisant I'egalitd : 
x = u 0 (7). 

[0024] Selon cette methode, il n"y a pas d'ambiguite a la mesure des vitesses parce qu'il n'existe qu'un seul maximum 
a cette fonction de correlation. 

so [0025] La methode de decalage temporel derive du deplacement r6el des cibles entre deux tirs successifs dans le 
temps selon le meme axe Z. C'est pourquoi la fenetre de resolution doit etre habilement choisie. La FIG.3A illustre 
une premiere fenetre de resolution de dimensbn w 0 = 2 u.S appliqu6e a un signal temporel d k (r,t) et la FIG.4B illustre 
une seconde fenetre de resolution de dimension w 1 = 4 u.s. appliquee au signal ulterieur d k+1 (r,t). La FIG.4C illustre la 
correlation C des signaux temporels d k (r,t) et d k+1 (r,t) qui conduit a la determination d'un pic C^xde correlation pour 

ss x de I'ordre de 0,3 u>s ce qui permet de calculer fa vitesse V z (r s t) par la relation (6). 

[0026] Afin de s'affranchir des erreurs inhSrentes a rechantillonnage dans la determination du maximum de la fonc- 
tion de correlation, le systeme de traitement 300 comprend favorablement un circuit 303 de multiplexage-interpolation 
fournissant a partir des valeurs des fonctions de correlation une estim6e plus precise de valeur du pic de correlation 
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et done de la Vitesse correspondante. Dans ce type de traitement, la correlation entre signaux est une correlation dite 
«1bit » en ce sens que les signaux y k+1 (r,t) et y k (r,t) utilises prec6demment sont reduits a leur signe. On sart que, dans 
ce cas, le pic C MAX de la fonction de correlation est de forme triangulaire isoc&le. Comme illustr§ par la FIG.4. la 
connaissance de cette forme permet de n'effectuer le calcul de correlation qu'en des points discrets des signaux tern- 
5 porels, et, une fois determines un certain nombre de points, permet a parlir du point le plus haut et de ses deux voisins 
de reconstituer completement, par interpolation Iin6aire, le picde correlation C MAX etdonc de determiner avec precision 
('emplacement de u 0 = t, d'ou Ton peut deduire V z (r,t). 

[0027] En reference a la FIG. 1 , les valeurs trouv^es pour la vitesse V(t,z) sont stockees dans une memoire 304 de 
stockage en vue de traitements ulterieurs. Comme illustre par la FIG. 5, en reportant les valeurs de vitesse V z (r,t) en 
10 cm/s en ordonnee et les profondeurs r en mm en abscisse, on obtient un profil des vitesses sur I'axe Z dans Partere 
au temps t. De meme que I'axe Z est discr6tis6 en N points d'une part de O et en N points d'autre part de O, le profil 
des vitesses est discretise avec le meme nombre de points. 

[0028] En reference a la FIG.1, a I'aide de moyens de visualisation connus, il est possible d'afficher sur un ecran 
305 une image en couleur de recoupment du sang qui est derivee de revolution au cours du temps du profil de vitesse 

is danslevaisseau. Lecodagecouleurde la vitesse dans cette image est noteeCVI (del'anglais ; Color Velocity Imaging). 
[0029] Un avantage important de la methode echographique de correlation en 1 bit est que les profils des vitesses 
sont obtenus en temps reel. Cependant toute autre methode apte a fournir des donnees de vitesse du sang dans un 
vaisseau avec une precision suffisante peut etre utilisee pour mener a bien revaluation de la viscosite et du gradient 
de pression comme decrit ci-aprfcs. 

20 [0030] Les profils de vitesse du sang, a la meme localisation Z, varient dans le temps. Les caracteristiques d'un 
fluide telles que la viscosite, la compressibility I'isotropie; la nature du flux telle que laminaire ou turbulente; les 
proprietes du tube qui contient le fluide, telles que la rigidite, I'elasticitS, ou la viscoelasticite ; et la forme du gradient 
de pression telle que large-bande ou bande etroite : sont des parametres qui dictent la forme des profils de vitesse. 
[0031] On suppose ci-apres que le sang dans uneartere est assimilable a un fluide Newtonien ayant un flux laminaire 

25 s'ecoulant dans un tube rigide. Comme illustrG par la FIG. 5, sous un gradient de pression note dPi dx, la vitesse V(r, 
t) suit une bi de Navier-Stokes, exprimee par la formule connue (8) : 



30 
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dP .dV. ,Gf*K 1 dV 



ou p est la densite du fluide et ou u, est la viscosrtd 

[0032] La loi de Navier-Stokes, qui lie la viscosite et le gradient de pression a la vitesse, ne permet cependant pas 
de calculer directement la viscosite du fluide et le gradient de pression en fonction de la vitesse du fluide. Le calcul de 

35 la viscosite du fluide par inversion de I'equation de Navier-Stokes conduit a des expressions analytiques complexes. 
[0033] Selon ('invention, ces difficultes sont evitees, et la viscosite et le gradient de pression sont determines en 
temps reel par un precede de traitement de donnees incluant une inversion de la relation (8) qui consiste a estimer la 
viscosite et le gradient de pression a partir des donnees fournies par des profils de vitesse prealablement determines 
de preference comme ddcrit pr6c6demment. 

40 [0034] Sous I'hypothese d'un flux laminaire, I'equation de Navier-Stokes gouveme les variations de la vitesse V(r,t) 
pour r et t donnes. Pour les points discr6tis6s Rj de R 0 a R N , de coordonn6es spatiales rj de r 0 a r N sur I'axe Z, 
I'integration de liquation (8) de Navier-Stokes entre r = r, et r = r n foumit I'expression (9) suivante, la d6rivee radiale 
de V x etant nulle : 



nrj dt " dx H /? dr K 1 



ou Tj est la profondeur d'exploration 
50 [0035] Chaque membre de I'equation (10) est une fonction du temps t et du rayon rj. Pour chaque echantillon de 
temps, la relation peut dtre ecrite sous une forme vectorielle, chaque vecteur 6tant forme par les valeurs des differentes 
quantites a differentes profondeurs d'exploration r y D'ou : 
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(10) 



L' equation du probleme direct peut done etre mise sous ia forme : 

y = u,h + 5 + b 



(11) 



ou b est un bruit blanc gaussien qui modelise les erreurs de mesure. En accord avec I'equation (11 ), y est un vecteur 
gaussien dont la moyenne est egale a : 



u.h + 5 



(12) 



[0036] La phase essentielle du procede consists a inverser cette equation qui modelise le probleme direct (11). La 
solution du probleme inverse est obtenue de preference au moyen d'un estimateur de maximum de vraisembiance 
25 connu (Maximum Likelihood Estimator) utilise pour estimer a la fois la viscosrte u. et le gradient de pression 5. 

[0037] On cherche les valeurs u. et 5 qui maximisent la vraisembiance de y connaissant u. et 8 donnees par une 
formute (13) : 



30 



(fij) - argmax (fl S) [p(y\/j,S)] (13) 



ou p(ylu., 6) est la probability de y connaissant n, 8, et ou u. et 8 sont des estimateurs de u. et 8 appelees estimateurs 
du maximum de vraisembiance. Ceci est equivalent a rechercher le minimum d'une fonction de cout notee W(j! ; 5): 

35 

W(M) = (y - uh - 8) T K" n (y -\ih - 5) (14) 

ou K" 1 n est une matrice diagonale dont les coefficients dependent de la confiance que Ton accorde aux valeurs de y, 
40 ce qui conduit aux expressions suivantes de |i et 5 : 



45 



y K n h — 

a = VK a J_ 

^ j -1 t -1 

h k n h — 

1 *«1 



(16a) 



50 



U h T l£y-Qh r dh {1?a) 



55 [0038] En premiere approximation, on suppose que le bruit est distribue de facon identique, de maniere a ce que la 
matrice K n soit proportionnelle a la matrice identite. Ceci produit une premiere approximation de la solution : 
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y h- — j— 



{1(0) = ^ (16b) 

s h h = — 

1 '1 

5(0)= ^-^ (17b) 
10 hi 

[0039] Selon invention , on propose un procede de determination de la viscosity u. fonde sur cette methode de 
resolution du probleme inverse, a partir du probleme direct constitue par la relation (11) et des profits de vitesse de- 
termines de preference par echographie, ou par toute autre methode donnant des resultats equivalents, en utilisant 

75 cet estimateur de maximum de vraisemb lance selon des etapes decrites ci-apres. 

[0040] En reference aux FIG. 1 et 6, un processeur 400 est assocte au systeme de traitement 300 de I'appareil 
echographique ou d'un systeme foumissant les vitesses. En echographie, le systeme de traitement 300 est «materiel 
» au sens HARDWARE utilise par I'homme du metier, et le processeur 400 est programme au sens SOFTWARE 
egalement connu. Le processeur 400 peut faire partie de I'echographe ou bien d'une station de travail couplee a 

20 I'echographe. La programmation du processeur 400 definit les differentes etapes n6cessaires pour determiner la vis- 
cosite et le gradient de pression 8 selon le precede de I'invention ilfustre par les blocs lonctionnels de la FIG.6. Ce 
precede comprend des etapes de : 

acquisition 410 des profils de vitesses V x (r,t) a partir des profils de vitesse V z (r,t) foumies par le systeme de 
25 traitement 300, en utilisant la formule (1 ) ; 

echantillonnage 420 d'un profil de vitesse discretise donne, en utilisant un nombre n de points a R n parmi les 
N points de discretisation R 1 a R N , tel que 2 < n < N ; les n points de segmentation definissent des rayons ^ a r n 
du tube assimile au vaisseau sanguin, notes de maniere generate ij; 

determination 430 d'un vecteur de vitesses constitue des vitesses aux points de segmentations R-j a R n , ces 
30 vitesses etant notees V x (r-,) a V x (rp) et d'une maniere generate V x (rj), 

evaluation 440 des debits du fluide dans le tube, e'est-a-dire du sang dans le vaisseau, dans les sections de tube, 
definies par les rayons rj de ^ a r n , ce qui fournit respect ivement les debits notes Q(rj) de Q(r-,) a Q(r n ) ; 
calcul 450 du vecteur y selon la formulation (10) en effectuant des calculs des quantites 



35 



40 



45 



50 



pour rj de ^ a r n ; 



Ces calculs (20a) comprennent des sous-etapes de : 
determination du d6bit Q(rp calcuie sebn la formulation suivante : 



Q(rj)-y2/r ri VV(r ( ) (20a), 



determination des derivees temporelles des debits Q(rj) selon : 



dQ(r.) 

pour^der, ar n (20b), 
55 et multiplication des derivees temporelles des debits par les quantites respectives avec rj de n, a r n : 
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p/nn (20c) 



5 



ou p est la masse volumique du fluide supposee environ egale a 1 ; 

calcul 460 du vecteur h selon la formulation (10) en eflectuant des calculs des quantites : 



2 W^'i) 



10 r j dr 

qui sont les derivees axiales des vitesses pour chaque point Rj de R A a R N definissant les rayons ij ; 
calcul 470 d'un produit matriciel : M T K" 1 M, ou la matrice M est donn£e par la formule : M= y - \ih - 8, et ou M T est 
la matrice transpose de M. 
15 calcul 480 de la viscosite u. selon les formulations (16) 

calcul 490 du gradient de pression 8 = -*fcselon (17) 

[0041] Pour effectuer les Stapes 410 a 490, le processeur est programme par les formules de derivees temporelles 
de debit, de derivees axiales de Vitesse, et de calcul matriciel. Le systeme de traitement 300 fournit automatiquement 
20 en temps reel les profils de vitesse et le processeur 400 fournit automatiquement subsequemment la viscosite et le 
gradient de pression dans un delai total extremement court c'est-a-dire en temps reel. 

[0042] Grace au procede selon I'invention, pour resoudre I'equation (8) qui a deux inconnues u. et 5, on procede" 
done par une inversion du probleme direct qui utilise n points discretises d'un profil de vitesse pour resoudre en realite 
n equations a deux inconnues, avec n > 2. Dans le cas de la determination de la viscosite du sang, il est interessant 
25 de selectionner un nombre n > 2 de points de discretisation pour disposer de plus de 2 equations independantes a 
deux inconnues pour la determination des 2 inconnues. En effet, si Ton se reduit a n = 2, on peut se trouver occasion- 
nellement devant une indetermination qui ne permet pas de determiner les deux inconnues. Dans ces conditions les 
moyens de calculs sont bloques. Mors que si Ton choisit favorablement un nombre n>2 de points de segmentation, il 
se trouve toujours un nombre suffisant d' equations pour determiner les deux inconnues u. et 5. 

30 



Revendications 

1. Procede de determination de la viscosite et du gradient de pression dans un vaisseau sanguin, comprenant des 
etapes de : 

acquisition d'une plurality n>2 de valeurs de vitesses du sang en correspondence avec un mfime nombre n 
de rayons du vaisseau sanguin, determines le long d'un diametre situe a une position axiale donnee, et a 
diffe rents instants, 

formation d'un vecteur de vitesses du sang au moyen de ladite pluralite de valeurs de vitesses du sang, 
evaluation de ladite viscosite et du gradient de pression a partir d'une transformation dudit vecteur de vitesses 
du sang. 

2. Procede selon la revendication 1 , comprenant, pour effectuer la transformation dudit vecteur de vitesses du sang 
des etapes de : 

formation d'un vecteur des debits (y) constitue d'une m§me pluralite n>2 de valeurs de debits du sang a travers 

des sections du vaisseau sanguin en correspondence avec les n valeurs de rayons, 

formation d'un vecteur (h) des d6riv6es des vitesses constituS d'un quotient de la meme pluralite n>2 de 

derivees radiales des vitesses par les rayons correspondents, ce vecteur etant factorisS par la viscosite (u,) a 

determiner, 

formation d'un vecteur de gradient de pression (8) constitu§ du produit du gradient de pression a determiner 
par une matrice identite, 

formation d'une relation lineaire liant directement le vecteur des debits (y) au vecteur des derivees de vitesses 
(h) factorise* par la viscosite (u.) et au vecteur de gradient de pression (8), 

evaluation simultan£e des deux valeurs a determiner de viscosite (u,) et de gradient de pression (8) a partir 
de cette equation directe. 
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Procede selon la revendication 2, comprenant revaluation simultanee des deux valeurs a determiner de viscosite 
(u.) et de gradient de pression (5) par une operation d'inversion de la relation lineaire directe liant le vecteur des 
debris (y) au vecteur des derivees de vitesse (h) factorise par la viscosite (u.) et au vecteur gradient de pression (8). 

Procede selon la revendication 3, comprenant, pour effectuer I'inversion de la relation lineaire directe liant le vecteur 
des debits (y) au vecteur des derivees de vitesse (h) factorise par la viscosite (u,) et au vecteur gradient de pression 
(5), une etape d'estimation du maximum de vraisemblance du vecteur des debits (y) connaissant la viscosite (ji) 
et le gradient de pression (5). 

Procede selon la revendication 4, comprenant, pour estimer ledit maximum de vraisemblance du vecteur des 
debits (y) connaissant la viscosite (u.) et le gradient de pression (5) des etapes de : 

evaluation a partir de ladite relation lineaire directe d'une premiere matrice (M) telle que: 

M = y - jin - 8 

et d'une seconds matrice transpose de la premiere matrice (M T ), 
evaluation d'une troisieme matrice (K) qui est une matrice diagonale dont les coefficients dependent de la 
distribution du bruit, 

et evaluation dudit maximum de vraisemblance du vecteur des debits (y) connaissant la viscosite (u) et le 
gradient de pression (8) en effectuant une recherche du minimum d'une fonction de cout (W) qui est le produit 
de la premiere, troisieme et seconde matrices tel que : 

MK 1 M T 

Procede selon la revendication 5, ou la troisieme matrice diagonale (K) est une matrice identrte diagonale, le bruit 
etant considere uniformement reparti. 

Procede selon la revendication 2, ou !e nombre des valeurs de vitesses est deux (n-2) et ou ('evaluation des deux 
valeurs a determiner de viscosite (u.) et gradient de pression (8) comprend des etapes de : 

formation d'un couple de relations lineaires directes liant le vecteur des debits (y) au vecteur des derivees 
de vitesse (h) factorise par la viscosite (u.) et au gradient de pression (5), de la forme 

y = uti + 8 

resolution du couple de relations lineaires directes, au moyen des deux valeurs de debits et de deux valeurs de 
derivees des vitesses, pour foumir la viscosite (u) et le gradient de pression (5). 

Procede selon Tune des revendications 1 a 7 comprenant pour I'acquisition de valeurs de vitesses des etapes de : 

acquisition de signauxtemporelsfoumis par un transducteur ultrasonore explorant un vaisseau sanguin selon 
un axe de tir echographique (2), 

traitement des signaux temporels dans le domaine temporel pour foumir des signaux correles temporellement 
et un profil de vitesses du sang discrelise le long de la direction de I'axe de tir, 

estimation d'un vecteur de vitesses du sang a partir d'un nombre n > 2 de valeurs de vitesses selectionnees 
sur le profil des vitesses a differentes profondeurs d'exploration (r) sur I'axe de tir (Z) dans le vaisseau sanguin. 

Procede selon la revendication 8, ou I'etepe de traitement des signaux echographiques temporels est une corre- 
lation temporelle en 1 bit. 

Appareil echographique pour la mesure et I'affichage de parametres physiologiques du sang s'ecoulant dans un 
vaisseau sanguin, incluant des moyens de mesure de la viscosite et du gradient de pression du sang, ces moyens 
comprenant : 

des moyens (200) pour I'acquisition de signaux temporels fournis par un transducteur (100) ultrasonore ex- 
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plorant un vaisseau sanguin selon un axe de tir echographique (Z), 

des moyens (300) de traitement des signaux temporels dans le domaine temporel pour fournir des signaux 
correles temporellement et un profil de vitesses discretisd le long de la direction de I'axe de tir (Z), 
des moyens (400) pour : 

5 ('estimation d'un vecteur de Vitesse du sang a partir d'un nombre n > 2 de valeurs de vitesses selectionnees 

sur le profil des vitesses a differentes prolondeurs d'exploration (r) sur I'axe de tir (Z) dans le vaisseau sanguin, 
et Evaluation de la viscosity du sang (u) et du gradient de pression (8), par une transformation dudit vecteur 
des vitesses en utilisant un proced6 selon Tune des revendications 2 & 9. 
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